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DISCUSSION 

of the Bates- Watts paper, 
Multiresponse Estimation With Special 
Applications to First Order Kinetics. 

Michael Frenklach 

Department of Materials Science and Engineering 
Pennsylvania State University l 



The authors presented an interesting approach to pa- 
rameter estimation for first order kinetic systems. The 
method is user oriented and particularly suited for com- 
puter implementation as a "canned" program. Indeed, 
present chemical kinetic codes input reaction mech- 
anism in a natural chemical language, that is, specifying 
reactions (usually in unformatted READ routines) as 
they are conventionally written on the paper. This infor- 
mation is automatically converted to a so-catted reac- 
tion matrix and, based on it, to differential equations 
describing the kinetics of reaction species. The reaction 
matrix, which contains all the stoichiometry of the sys- 
tem, can conveniently provide the required input infor- 



1 Michael Frenkiach's contribution to the subject stems from work 
performed in the Department of Chemical Engineering, Louisiana 

State Univmtiy. 



Another important feature, from the user's point of 
view, is that the presented method is applicable to mul- 
tiresponse data. It should be realized that modern prob- 
lems of interest to chemical kinetics get tougher, as for 
example, formation of pollutants in hydrocarbon com- 
bustion. The experimental answer to the growing com- 
plexity of the systems is the employment of multiple 
diagnostics For simultaneous monitoring of various pro- 
cess variables. However, interpretation of the experi- 
mental results cannot be fully realized without reliable 
and convenient multiresponse methods. 

The following are some of my thoughts on the needs 
in this area: 

1) Often, kineticists exhibit a philosophical re- 
sistance to a multiparameter approach to experi- 
malion for automatic coding of the method of Bates and 
Watts. 
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mentation. A classical way is to "isolate" a given reac- 
tion of interest; under such conditions the rate 
coefficient parameters can be determined by a simple 
well-established straight-line treatment. Determination 
of more than one rate coefficient in a single set of experi- 
ments is considered "not clean experimentation." In 
principal, however, the isolation is not possible: there 
are always other reactions occurring simultaneously 
with the one of interest. The researchers usually engage 
in an elaborate line of reasoning to assume, sometimes 
unjustifiably, single-reaction conditions. These kinet- 
icists must realize that multiparameter analysis using rig- 
orous multiresponse techniques can provide more accu- 
rate and informative answers. Neglecting, for instance, a 
chemical reaction with the rate contribution of, let us 
say, 10%, can lead to a much larger than 10% distortion 
in the estimation of the main parameter. Statisticians, on 
the other hand, should demonstrate the techniques they 
develop on examples of current interest and difficulty. 

2) Although first order kinetic models constitute an 
important class, higher order kinetics are of more gen- 
eral interest and there is a great need for development of 
statistical methods for these nonlinear systems. 

3) Most estimation methods, including the one 
presented by Bates and Watts, concentrate on deter- 
mining the solution which minimizes the objective func- 
tion and only approximate confidence limits. What is of 
interest to many applications is the joint confidence re- 
gion. It should be noted that in the problems of chemical 
kinetics these regions are usually not ellipsoidal, for 
which second order approximation methods are suf- 
ficient, but crescent shaped. 

4) While estimating parameters, it is most important 
to check the model adequacy. This point was excellently 
demonstrated by Box and Draper (1965). These authors 
warn that "the investigator should not resort immedi- 
ately to the joint analysis of responses. Rather he 
should... consider the consistency of the information 
from various responses." To my knowledge, however, a 
formal multivariate lack-of-fit test for a general non- 
linear case has not been developed. 

5) A question on the number of degrees of freedom 



was brought up by Bates and Watts. Using fast 
digital sampling electronics, the number of observations 
per response can be very large (in our laboratory this 
number was approximately 1000). Does this number de- 
termine the degrees of freedom? If so, then one can 
easily increase this number by orders of magnitude by 
using faster electronics. This point should be clarified. 
Finally, I would like to point out that in an attempt to 
resolve some of the issues brought up above, a method 
for multiresponse parameter estimation applicable to a 
dynamic model of general order was developed in our 
laboratory (Miller and Frenklach, 1983; Frenklach, 
1984; Frenklach and Miller, 1985). The method is based 
on approximating the solution of the differential equa- 
tions describing the kinetics of reactive system instead of 
the equations themselves. The approximation is devel- 
oped following the methods of empirical model building 
(Box et al. 1978) and the concept of computer experi- 
ment of Box and Coutie (1956). Once the approxi- 
mations to all responses are obtained, the parameter esti- 
mation, determination of joint confidence region, and 
lack-of-fit test are easily performed following the ap- 
proach of Box and Draper (1965). 
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